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q\ ■ Abstract 

Oh! 

■ , It would be interesting to investigate the accuracy of the results obtained 

in the variational method, because it is important for studying hadron spec- 
if ' tra. One can define some criteria to judge the accuracy, or the quality of the 
trial function. We employ a simple potential form to check how accurate the 
S^ . variational results obtained by a single-parameter trial function can be. All 
the concerned problems, in particular, the relevant aspects on the applica- 
tion of hypervirial theorem in variational method for various potential forms, 
are discussed in every detail. The qualitative conclusion of the work can be 
generalized to much more complicated cases. Our study suggests that the 
hypervirial relations can serve as a good and practical criterion for judging 
the accuracy of any trial functions. 

PACS number(s): 03.65.Ge, 11.80.Fv, 12.39.Jh, 12.39.Pn 
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I. INTRODUCTION 

As well known, there are few potential forms whose bound state problems can analyt- 
ically be solved. Therefore, various approximate methods have been used to do the job, 
in which, no doubt, numerically solving the Schrodinger equation is the most direct and 
powerful method, which in principle can meet any desired accuracy. But on other side, 
it cannot provide analytical expressions for further discussions. Moreover, usually, as one 
numerically solves the Schrodinger equation, he can only keep a few terms in the zero-th 
order, while leaving the others, such as the spin- dependent terms to be treated as a pertur- 
bation. Indeed, the perturbative method is an approximate treatment which has extensively 
been adopted because of its simplicity and convenience. However, this method is limited by 
the convergence of perturbative expansion. In fact, treating wave functions is much more 
difficult than evaluating energy eigenvalues [0. 

Another widely employed approach is the variational method. Most desirable advantage 
of the method is that there exists an analytical expression of the wave function for further 
discussions on physics. In past years, lots of work paid their attention to seeking for accurate 
energy eigenvalues but seldom discussed the wave function. In many simple cases, the 
Gaussian-type or the exponential-type trial wave functions (TWF) with a single variational 
parameter are employed and they usually can give an accurate value for the eigenenergy, 
but it is not true for the wave function in general. 

In a recent work [F2], we investigated this problem in detail with various potential forms 
which are used frequently in literatures to study heavy quarkonia and put more emphases 
on estimating the wave function at origin (WFO) by using a TWF with single and several 
variational parameters. We hope that it would shed some light on how to obtain sufficiently 
accurate WFO in terms of the variational method. We find that the accuracy of WFO 
seriously depends on the choice of the forms of TWF. In general, although there is no 
universal rule to determine the form and the number of the parameters of the TWF, one 
can always find out a relatively simple and more reliable TWF with the least number of 
parameters. Our result showed an economic way to realize it. 

Recently, Lucha and Schoberl have made a renewed investigation || to qualitatively 
and quantitatively evaluate the accuracy, namely how close to the exact eigenstates the 
results obtained in variational method can be. They embark upon a systematic study of 
the accuracy of the variationally determined eigenstates of a Hamiltonian H and suitable 
measures to judge their quality and apply their principles to the prototype of all (semi- 
) relativistic bound-state equations, the "spinless Salpeter equation". In their work, the 
harmonic-oscillator potential is employed as an example. Their work is very instructive and 
stimulates interest in studying the accuracy of the trial wave function for the variational 
method. It is the goal of this work to further study the accuracy of the TWF with single 
variational parameter and the criteria to judge the accuracy. 

In this work, we firstly extend the criteria used by the authors of Ref. || and apply 
them to treat a simple case of the non-relativistic Schrodinger equation where only the 
linear potential exists and give some further discussions on the accuracy of the TWF with 
single parameter. The motivation to choose the linear potential is three-fold. The first is 
its significance in the particle physics. The next is that if only the linear potential exists, 
there are exact solutions for nS states, which are the well-known Airy function Ai(z). Thus, 



we can conveniently use them as a basis for comparing the results obtained with the TWF 
and the exact one to quantify how close to each other the approximate solutions in the 
variational method and exact eigenstates can be. The last one is that we find that the linear 
potential possesses some interesting features which are helpful for a thorough investigation 
of the criteria which can be generalized to other cases. 

After the quantitative discussion to all the criteria for the ground states and the lowest 
radial excitation state of the linear potential, we study in detail the application of hypervirial 
relations as a most powerful and practical criterion to judge the accuracy of the variational 
trial wave function. We shall prove that when the variational parameter a take its optimal 
value, the special form of the TWF 

R&W) = N ^ e ~ a rb 

satisfies the general virial theorem which is so called the first order hypervirial relation. 
This is independent of both the concrete form of the central potential and the specific value 
of b. In this case we have to use higher order hypervirial relations as a criterion to judge 
the quality of the TWF. To further testify the applications of the hypervirial theorem, we 
calculate the relevant quantities corresponding to several commonly adopted potential forms 
such as the Cornell, logarithmic, Martin's and power-like ones. 

The paper is organized as following. After the introduction, we provide several criteria 
for judging the quality of the TWF, where the linear potential is taken as an example to 
serve as the basis for discussion. Indeed applications of all of them need to priori know 
the exact solutions for a comparison, studying them is to find support to our main goal, 
i.e. to apply the hypervirial relations as criteria. In Sec. Ill, we discuss applications of the 
hypervirial relations which can serve as practical and useful criteria for the quality of the 
TWF. In the last section, discussion and conclusion are made. 



II. SEVERAL EXACT-SOLUTION-DEPENDENT CRITERIA FOR TWF 



What we concern is how close to the accurate solution the obtained approximate wave 
function can be when the approximate energy E m i n of Eq.([B5| in the Appendix B) is sat- 
isfactorily close to the accurate energy eigenvalue, or in other words, how seriously the 
approximate Rtwf deviate from the real one and how the nature and extent of these errors 
affect the calculated values. There were intensive discussions on the problem many years 
ago ||. Recently Lucha and Schoberl || reviewed and developed these discussions and 
introduced several measures to judge the quality of an approximate wave function in the 
variational method, whose applications to our case are described in the following. 

(a) The distance between the obtained approximate value of the energy E min and exact 
eigenvalue E. It can be represented by using the relative error e of E min and E, 
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E 
(b) The overlap integral S of the trial state \Rt w f) an d the eigenstate \R e i g ). 



(1) 



If both \Rtwf) and \R e i g ) are normalized, i.e. 

\Rtwf\Rtwf) — \Reig\Reig) = 1) 

one can define 

S = (Rtwf\Reig) ■ (2) 

If the R tw f is just the exact eigen-state R e i g , the overlap S is equal to unity. In general, we 
have 

0< \S\ < 1. 

As a quantitative criterion, according to the suggestion of authors of |§, we take the deviation 
of the squared modules of the overlap 5*, from unity, a as 

a = l-\S\ 2 . (3) 

(c) The normalized maximum difference of variational function R tw f(r) and the eigenstate 
R e i g (r), i.e., the maximum relative error, 

max[R twf (r) - R eig {r)} 

u = — — . (4) 

max R eig {r) 

Besides the above criteria, we introduce the following measures which may reveal the 
physical properties of the TWF. 

(d) The squared wave function at origin (WFO) R 2 (0). 

As is noticed in our earlier work |2| , this quantity can be obtained by using two methods. 
The first approximate WFO, R 2 (0) I, is directly calculated from the normalized radial wave 
function by setting r = 0. And in the second approach, R(0) 2 II can be obtained from the 
well-known equation || 

i? 2 (0) = 2/i( R\— \R). (5) 

If R(r) is the real solution of the Schrodinger equation with the given potential, the resultant 
-R(O) 2 I and R(0) 2 II should exactly be the same. However, if we use R tw f(r) to calculate 
these two quantities, the results definitely deviate from each other. Our numerical results 
H showed that the later one can generally reach very high accuracy. 

In the present case of the linear potential V(r) = r with 2/i = lGeV, it is very easy to 
prove that the resultant R(0) 2 II is 

R 2 (0) = (R\R). (6) 

Thus, for the any R(r), no matter it is an accurate eigen-wavefunction R eig {r) or any TWF 
Rtwf( r ), so long as it satisfies the normalization condition, the i?(0) 2 obtained by using such 
method must be identical to unity. Definitely, this is a special feature possessed by the 
S'-wave bound state for the linear potential. Therefore, we only have the first type of the 
WFO R(0) 2 I can be used in our case. In practice, what we take as a quantitative criterion 
is the relative deviation 



r R twf{fy ~ R lig{fy /~s 

^ S 76,(0) ■ (7) 

However, in the general case, dV/dr ^ 1, then both i? 2 (0)I and i? 2 (0)II can be used to 
obtain 5 ro 2- 

(e) The average value of the operators r 2 and r~ 2 . 

In order to show the effect of the approximate wave function on some particular prop- 
erties, we calculate average values (r 2 ) and (r~ 2 ). In fact, the weighted overlap integrals 
H depend principally on regions of the configuration space. The accuracy of (r 2 ) depends 
mainly on the accuracy of TWF in the outer regions and that of (r~ 2 ) depends more sen- 
sitively on the region near the origin. In practice, what we select as the criteria are the 
relative errors of the approximate value with respect to the exact one 

r _ (Rtwf\ r \Rtwf) ~ (Reig\ r \Reig) , Q s 

K2 = (R~F\R~) (8) 

and 

f \Rtwf\T~ \Rtwf) — (Reig\r~ \R e ig) /riX 

6rn2 ~ (R ■ \r- 2 \R ) • (9) 

\ J - L eig\' \ J - L etg/ 

All the numerical results corresponding to these criteria for the ground state and low 
excitations of the linear potential, including the values of the optimal variational parameter 
ao, the minimum values of the energy E m i n , are listed in Table 1. 

The numerical results in Table 1 confirm that all the criteria can be effectively used for 
judging the quality of the TWF. 



III. THE HYPERVIRIAL RELATIONS AS MOST PRACTICAL AND 

POWERFUL CRITERIA 

A. The hypervirial theorem and hypervirial relations 

The virial theorem has been widely used in both classical and quantum mechanics. In 
quantum mechanics, it can be formulated as 

2 <T >=<f • W(r) >, 

where T is the kinetic energy operator and the expectation value is taken for eigen- 
wavefunctions of the Hamiltonian H . It holds for any central potential V(r). 

It is easy to prove that an equivalent statement of the virial theorem is that the expec- 
tation value of the commutator [W, H] on the energy eigenstates \R e i g ), vanishes: 

{R eig \[W,H]\R eig ) = 0. (10) 

where 



W = i(r.p + p.r), (11) 

is called as the generator of dilatations and H is the Hamiltonian. 

Hirschfelder generalized the theorem to the so called hypervirial theorem 0. It states 
that the expectation value of the commutator [G,H], where G stands for any time- 
independent linear operator, on the energy eigenstate \R e i g ) is equal to zero: 

(R eig \[G,H]\R eig ) = . (12) 

Obviously our W is such an operator. 

However, for a trial state \R tw f), in general, the expectation value 
{Rt w f\[G, H]\R tw f) is not equal to zero. We can expect that the closer to the real solution 
the approximate trial function is, the closer to zero the expectation value would be. Hence, a 
set of operators G generates a whole set of the "hypervirial relations", which will be defined 
in Eq.(|TP|), and each of them may serve as a measure of the quality of a given trial function 



\Rtwf) by evaluating the expectation value {R tw f\[G, H]\R tw f). 

In the present work, the hypervirial operator can be chosen as J7LQ 

G = f(r)p r (13) 

where f(r) is a function of r and p r is the radial momentum operator conjugate to r, whose 
explicit expression is p r = —iidjdr + 1/r). The corresponding Hamiltonian (with 2/i = 1 
GeV) reads 

,2 , Kl + l) 



H=(PrY + ^^ + V(r). (14) 



Considering 



[G,H] = [f(r)p r ,H] 



and 



df(r) 2 d 2 f(r) ,, JU + 1) ,, ,dV 



* = B-*±*L-V (16) 



<?/(r) id?f(r) df(r) 

P' = ~— T3r-l—r-,H), (17) 



<ir 2 2 dr 3 dr 

we can get the commutator 



dr r 3 r 2 dr 2 <ir 3 

iy i d/(r 

dr 2 dr 



«/<r)?-i[^,fli. («) 



The hypervirial theorem demands the expectation values of both [G, H] and [-^p-, H] for 
the eigenstate of H to be zero, thus from Eq. (|T8| ) we can obtain an algebraic equation 



-<fWfW)£>-^ +1 )(<^>-<^>) 

where (■ • •) stands for the expectation value over the eigenstate. These algebraic equations 



(|i~9|) are defined as the hypervirial relations. 
Taking 

f(r)=r N (iV= 1,2,3,.-.), 

Eqs. C |T~8D and fll9|) can respectively be re-formed as 

[G N ,H] = [r N p r ,H] 

= 2iNr N -\H -V)- %r N< ^- - 2il(l + 1)(N - l)r N ~ 3 + -N(N - 1)(N - 2)r N ' 3 

dr 2 

- l -[Nr N -\H] (20) 

and 

2EN(r N - 1 ) = 2N(r N - 1 V) + (r N (DV)) + 21(1 + 1)(JV - l)^" 3 ) 

-l N (N-l)(N-2)(r N - 3 ), (21) 

where Gn = r N p r is the adopted hypervirial operators and Eq. (|2l|) is an important recur- 
rence relation which has extensive applications. 

For the nS state in the linear potential case, we can obtain, from Eq . ([19|) , a set of very 
useful recurrence relations 

2EN(r N ~ 1 ) = (2N + l)(r N ) - -N(N - l)(N - 2)(r N - 3 ), (22) 

which can be named the iV-th order hypervirial relation. 
When N — 1, we have the first order hypervirial relation 

2E = 3(r), (23) 

which is just the general virial theorem in the linear potential case. Taking N = 2 and 
iV = 3, we can obtain the second and the third order hypervirial relation, which respectively 
are 

AE(r) = 5(r 2 ) (24) 

and 

6E(r) = 7(r 3 ) -3. (25) 

In the above equations, the average is taken over the eigenstate of H . 



B. Deviation from the hypervirial relations as criteria for accuracy of TWF 

As mentioned in the last subsection, for a properly chosen trial function R tw f(r), these 
hypervirial relations should merely approximately hold. According to the 1-st, 2-nd and 
3-rd hypervirial relations, we can respectively define deviations: 

u 1 =2E-3{r), (26) 

z/ 2 =4£(r)-5(r 2 ), (27) 

and 

u 3 = -QE(r 2 } + 7(r 3 ) - 3, (28) 

where the average value of r N can be analytically calculated in terms of our TWFQ. For the 
ground state, by using Eq. ( |BT|) with the optimal variational parameter ao, one has 

r(M) (1 + 6) r(i) f 
(r > ~ r(§) I. 2r(|) > ■ (29) 

In terms of Eq.(|2P|) and E min , we can easily obtain the numerical values of U\, v<i and z/ 3 . In 
order to evaluate hypervirial relations for judging the quality of TWF, we further introduce 
relative errors 8vi (i=l,2,3...) instead of z/j as 

x - N x - I 272 ' x - M (?n\ 



Obviously, these expressions (|30|) of 5vi are only valid for the linear potential, for other 
potential forms, the denominators in (|30D should be replaced. 

The results of u±, 5vi, z/ 2 , 5z/ 2 , v$ and bv 3 are listed in Table 1, respectively. In the next 
section, we will give more discussions on the implications of these results. 

In this work, we specially choose the most frequently used hypervirial operators Gn = 
r N p r . When N = 1, G\ = rp r has a clearly physical meaning. Although it itself is not a 
Hermitian operator, it only differs from its symmetrized and Hermitian form W (see Eq . ([TT|) ) 
by a constant, i.e., 

W = G 1 - % -. (31) 

As a consequence, 

[W,H] = [G 1 ,H], 

and one obtains the same first order hypervirial relation. 



x It is noted that in this subsection, the <> stands for an average with respect to the trial 
wavefunction TWF 



As well known ||, the dilation generator W generates the scaling transformation or 
dilation of the phase space, which is an important transformation not only to quantum 
mechanics but also to classical mechanics. Its classical equivalent r • p has been used to 
derive a virial theorem in the classical mechanics. 

Using the hypervirial relation of the quantum mechanics derived in terms of the operators 
Gn-i one can derive many useful recurrence relations between (r N ) and E. These relations 
can be used to gain approximate eigen-spectra and the expectation values of observables 
without solving the wave functions of bound states as well as to calculate the phase shifts 
and the amplitudes of the scattering processes more effectively than in the perturbative 
method. In this investigation, what we concern is another important application noticed 
by Lucha and Schoberl ||. We will show that the hypervirial relations could be sensitive 
criteria for the accuracy of the solution by using the variational method. Of course, if we 
choose a more complicated form, we cannot have such compact expressions and the results 
of Pi and 8vi would be presented only numerically. Even so, the principle and the general 
conclusion hold. 

In the last section, we give many criteria which can be used to judge the accuracy of 
TWF. It can be seen from Table 1 that all these criteria give a consistent result, by which we 
can select a most optimal TWF corresponding to b — | for IS state of the linear potential 
from four possible value of b. However, the practical significance of these criteria to judge the 
accuracy of the TWF are limited by their dependence on the exact solution. It is clear that if 
one can get the exact solution, it, in general, would be unnecessary to find the approximate 
solution. But sometimes we want to use the variational method to find a simple analytical 
form of the approximate wave function with single parameter or a few parameters to replace 
a complicate exact solution or a pure numerical solution, and then these criteria are certainly 
very helpful. However, for those problems that we do not need an exact solution or it is 
very difficult to obtain the exact solution, these exact-solution-dependent quantities have 
no practical usage. In the case, as criteria the deviations from the hypervirial relations are 
the unique choice because only they do not require any knowledge on the exact eigen state 
and consequently can be applied to all the cases. 

Table 1 

The values of various criterion quantities for the single-parameter variational solutions of the 15*, 2S 
and 35* states in the linear potential case. The exact eigenvalues of energy are Eis = 2.33810Gey, E2S — 
4.08794Gey and E 3S = 5.52055GeV, respectively. 



Quantity 


IS state 


2S state 


35 state 


b 


1 


3 
2 


7 
4 


2 


7 
4 


7 
4 


a 


0.908560 


0.471405 


0.348957 


0.260531 


0.348957 


0.348957 


E m in (GeV) 


2.47645 


2.34723 


2.33825 


2.34478 


4.09917 


5.52874 


e 


0.05917 


0.00390 


0.00006 


0.00285 


0.00275 


0.00143 


a 


0.032731 


0.00207 


0.00003 


0.00180 


0.00634 


0.22584 


U) 


0.14303 


0.04353 


-0.00436 


0.35934 


-0.0808 


0.63515 


Sr 02 


2.00000 


0.33333 


0.03646 


-0.15117 


-0.12956 


-0.73780 


8r p2 


0.24648 


0.03070 


-0.00086 


-0.01264 


-0.09030 


-0.16022 


8r n2 


0.46777 


0.11295 


-0.00790 


-0.07351 


0.01229 


-0.30906 


V\ 














0.13701 


0.23532 


8v\ 














0.05278 


0.06820 


V2 


-1.81712 


-0.33349 


0.01439 


0.26760 


2.02267 


8.04999 


5v 2 


0.10000 


0.02220 


0.00099 


0.01859 


0.04989 


0.11795 


"3 


13.000 


1.9280 


-0.13021 


-1.5000 


-9.7304 


-53.865 


5u 3 


0.18571 


0.04000 


0.00298 


0.03571 


0.50498 


0.13402 



C. Restriction on applying hypervirial relations 

As noticed, there are some cases in which the expectation value (R tw f\[G, H]\R tw f) van- 
ishes accidentally. We list following possible cases: 

(i) TWF is an eigenstate of G, 

(ii) TWF is an eigenstate of the commutator [G, H] with eigenvalue 0, 

(hi) G\R tw f) is orthogonal to H\R tw f), 

(iv) \G,H]\R twf ) is orthogonal to \R tw f)- 

Evidently, in these cases the corresponding hypervirial relation cannot be used as a 
criterion. Table 1 indicates that we just confront this case for the IS* state in the linear 
potential. The optimal TWF ( [F5T| ) gives V\ = and the result is independent of the value 
of b. 

In fact, since we choose G\ = rp r as the first order hypervirial operator which is just the 
generator of the dilation W and has no Hilbert-space eigenvector ||, with N=l and Eq. (p0|) , 
we have 



[rp r ,H] = 3r-2H. 



(32) 



One can straightforwardly prove that for the optimal TWF of the ground state, R tw f(r) in 
the form of Eq. (plf), [rD, H]\R tw f) is orthogonal to the state \Rt w f)- Thus, this is the 
case (iv) and the 1-st order hypervirial relation cannot be used to judge the quality of the 
concerned TWF. Definitely, in this case, we can use v 2 which is not equal to zero and can 
sensitively judge the accuracy of the concerned TWF. To our knowledge, the authors of most 
literatures only discussed the possibility for the case where the hypervirial relation fails to 
testify the given trial function, but seldom gave concrete examples. The form of TWF ( |E>T| ) 
presented by us is a convincing illustration. 



10 



We find that at this point we can give stronger arguments and generalize it to extensive 
cases. Epstein and Hirschfelder have proved a well-known lemma [nj which states that if a 
variational trial function ip obeys 

Ta = M * (33) 

where a is the variational parameter and W is a Hermitian operator, then the optimal i/j 
satisfies the hypervirial theorem for the Hermitian operator W, i.e. 

(nw,HM= o. 



The authors of ref. [11] generalized the lemma and claimed that if Eq. (R3) is replaced by 



Ji = BiWip + iCip, (34) 

da 



the lemma is still correct. In Eq.(|34|) B and C are the functions of a , B is real and W is 
Hermitian. Evidently, Eq.Q3"3]) is a special case of Eq.(^) with B = 1 and C = 0. This 
lemma and its generalized form have extensively been used to select approximate TWF's 
and to discuss the spectroscopy of molecule and atom. 

Later, as a sample of the application of the generalized lemma, we will stay in the linear 
potential case, unless it is specially specified. 

With this generalized lemma, we can deduce some very useful consequences. 



(a) For the IS state, when TWF is in the form (|B 1|) , we have 



dR twf (r) 3 6 



da 2ab 

As aforementioned, taking 



7T- u -r b )Rt wf (r). (35) 



1 A Q 

W = -(rp r + p r r) = -i( r -T- + 7j) (36) 



as the hypervirial operator, we get 



iWR twf (r) = (~ - abr b )R twf (r). (37) 



Comparing ( 37|) with (35f), we can obtain 



dR twf (r) _ I 
da ab 



-iWR twf (r), (38) 



which corresponds to B = — -k and C = in Eq.(^). When the variational parameter a 
takes its optimal value, we must have v\ = 0, which is independent of the value b. This means 
that the optimal TWF ( |B1| ) for the IS* state satisfies the general virial theorem although it 
is not the exact wave function. If the employed TWF is of this kind, we at least have to use 
the second order hypervirial relation v 2 - 

(b) In fact, (a) can be generalized to the other cases, because the result is independent of 
the concrete form of the central potential V(r). As long as we take Eq. (|B1~| ) as TWF for the 

11 



IS state in any central potential case, the general virial theorem must be satisfied, namely 
v\ = 0. The general virial theorem cannot be taken as a criterion to judge the quality of 
TWF. Then, we have to use the deviations from 2-nd and 3-rd order hypervirial relations 
v 2 and z/ 3 as the criteria. 

To investigate the applications of the above consequence, we evaluate z/^s with the fre- 
quently used potential forms including (i)the Coulomb potential, (ii) the Cornell poten- 
tial V(r) = — + r, (iii) the Martin potential V(r) = r ai , (iv) the logarithmic potential 
V(r) = log(r) and (v) the three dimensional isotropic oscillator potential V(r) = r 2 , respec- 
tively. The expressions of z/j and <5z/j ( i=l,2,3) obtained from the general formula and the 
corresponding numerical results for these potentials are shown in Table 2. 

It can be seen from the table that for the cases , (ii) , (iii) and (iv), when b = |, one 
obtains the best TWF. On the other hand, in the case (i) (or (v)), the closer to 1 (or 2) the 
value of b is, the smaller the corresponding 5i>2 and Su^ are. Eventually, when 6=1 (or 2), 
the TWF is just the exact solution for the corresponding potential and all z/j and <5z/j go to 
zero. Obviously, one can use <5z/ 2 and 5v% as the criteria to judge whether TWF is accurate 
enough. 

Table 2 



Variational results with the single-parameter TWF (Bl) in some central potentials cases mentioned in 
the text. 
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b 


1 


3 

2 


7 
4 


2 


Coulomb potential V(r) = — - 


a 


0.50000 


0.19159 


0.11770 


0.07074 


E 


-0.25000 


-0.23555 


-0.22400 


-0.21221 


z/i = -(1/r) -2£ 














i/ 2 = -3 - AE(r) 





-0.31289 


-0.39463 


-0.45352 


i/ 3 = -5(r) -6£(r 2 ) -3 





-3.1523 


-3.9158 


-4.5000 


^ 2 -¥ 





0.10439 


0.13154 


0.15117 


<^= a 





0.01543 


0.01931 


0.02122 


Cornell potential V(r) = —1/r + r 


a 


1.10939 


0.63544 


0.49016 


0.37969 


£ 


1.47345 


1.39969 


1.41824 


1.45064 


i/x = 3(r) - (A) - 2E 














v 2 = 5(r 2 ) - AE{r) - 3 


1.21878 


-0.08892 


-0.40439 


-0.63714 


1/3 = 7(r 3 ) - 6£(r 2 ) - 5(r) - 3 


7.1409 


-0.35635 


-1.8015 


-2.7949 


S"2 - g (r 2) 


0.10000 


0.00881 


0.04094 


0.06514 


<^3 - J^ 


0.18571 


0.01370 


0.07375 


0.11708 


Martin potential V(r) = r ai 


a 


0.24300 


0.06770 


0.03641 


0.01965 


£ 


1.24007 


1.23576 


1.23654 


1.23811 


1/1 = 2.1(r ai ) -2£ 














z/ 2 ^^(r 1 - 1 ) -4E(r) 


0.26731 


-0.01636 


-0.08947 


-0.14441 


1/3 = 6.1(r 3 ) -6E(r 2 ) -3 


6.6550 


-0.41743 


-1.8517 


-2.8500 


x,, — 1^1 


0.00865 


0.00058 


0.00320 


0.00515 


^2 - 4.^1.1) 


X ; , — 1^3 I 


0.01717 


0.00140 


0.00645 


0.01005 


0Z ^ - 6 . 1(r ,3) 


Logarithmic potential V(r) = log(r) 


a 


0.70711 


0.33693 


0.23663 


0.16667 


E 


1.07521 


1.0450 


1.05310 


1.06755 


v x = 2{log(r)) + 1-2E 














v 2 = 4(r %(r)) - (4E - l)(r) 


0.70711 


-0.08362 


-0.28942 


-0.44455 


1/3 = Q(r 2 logr) - (QE - l)(r 2 ) - 3 


6.0000 


-0.67223 


-2.0412 


-3.0000 


Si' — ^ 2| 


0.09162 


0.01362 


0.04854 


0.07474 


ui -2 — 4 ( r /o9 ( r )) 


5i/ — '^ 3 


0.14374 


0.02480 


0.08129 


0.12334 


""■* 6<r 2 logr> 


3-D isotropic oscillator potential V(r) = r 2 


a 


1.3161 


0.78082 


0.52127 


0.50000 


£ 


3.4641 


3.0667 


3.0136 


3.0000 


i/! = (r 2 ) - E/2 














v 2 = 6(r 3 ) -AE(r) 


3.9482 


1.0499 


0.43243 





u 3 = 8(r 4 ) - 6£(r 2 ) - 3 


21.000 


4,63819 


1.81366 





Su 2 - 6{r 3) 


0.37992 


0.11412 


0.04783 





5^3 - 8(r4) 


0.35000 


0.12937 


0.05657 
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D. Application to the IP case 

Now, let us turn to the IP state. The wavefunction for the IP state is Ri p (r)Yi r 
For the linear potential case, the differential equation for R\ p {r) reads 

, d 2 2d, 2 



1 dr 2 r dr r 2 



R lp (r) = E R lp (r). (39) 



There is no analytical solution for Ri p (r). Although one can solve it numerically, if we use 
the variational method to find an optimal TWF for the IP state, we would not have a 
convenient exact solutions like Airy function for the S state to compare with. In this case, 
only the hypervirial relations can serve as a powerful measure to judge the quality of TWF. 
Now, we write the optimal TWF in a more general form 

I$2f{r) = N n r n e—* (n = l), (40) 

with the normalization constant N n being 

, . 3+2n , . , 

N n = {2a)^r .f^j (n = 1). (41) 

It is very easy to prove that this TWF also satisfies the general virial theorem, namely 
fi = 0. Because we have 



ggS(r) _ 3_+2n 6 (n) 

da ~ { 2ab r > K ^ r > ^ 



and 



comparing (^) with ((4~3"|), we can immediately reach Eq. fl38|) . Actually, R tw f(r) is a special 
case of TWF R$ f {r) with n = 0. 

Taking n — 1, one has the IP state TWF 

R W f{r) = Nire ~ar \ (44) 

We can fix the optimal values of the variational parameter a and the energy E for b = 1,1,1 
and 2, respectively, according to the normal procedure, and the results are shown in Table 
3. The analytical expressions and their numerical results of the deviation from the first and 
second hypervirial relations U\, v-i and 5v2 are also shown in the table. 
It is easy to see from Table 3 that 

(a) v\ is indeed zero, which is independent of the value of b. This is consistent with our 
general discussion on the TWF P^/(r); 

(b) The optimal value of a is different from the one for IS state, and 8v 2 can sensitively 
be used to select the best TWF. The result is b = ? again. This conclusion is consistent 
with that achieved for the IS state. 

Table 3 

Variational results with the single-parameter TWF Eq.(|4J) for the IP state in the linear potential case. 
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b 


1 


3 
2 


7 
4 


2 


a 


1.07722 


0.47141 


0.32323 


0.22452 


E 


3.48119 


3.36927 


3.36135 


3.36778 


v x = 3(r) - 2P 














V2 


= 5(r)+4(i) -AE(r) 


2.15443 


0.45059 


-0.02146 


-0.39152 




St — ^ 2 


0.06250 


0.01467 


0.00071 


0.01311 


dl/ * ~ 5<^>+4(I) 



Instead, if we still employ a TWF in the form of 



Ri P {r) 



Core 



7 



(45) 



where a$ is taken to be the optimal value for the IS state and Co is determined by the nor- 
malization condition. Then, there is no free parameter. Calculating (H), one can straight- 
forwardly obtain the energy 

E = 3.6788. 

The corresponding deviations from the first and second order hypervirial relations are 

Vl = -0.1003 Sui = 0.04676, 
and 



V2 



-1.0472 5u 2 = 0.03760, 



respectively. This is because that the Hamiltonian for the IP state is different from that 
for the IS state, so that the adopted ao is not the optimal value for the IP state, and the 



Pip(r) given by (45) does not satisfy the general virial theorem, namely z/i is not equal to 
zero. 

In comparison with the values of U2 and 5^2 in Table 3 for b — |, one can conclude that 
the TWF chosen in such a way is not an appropriate one. 

Now, we add one more variational parameter in and take TWF in the form of 



Ri p (r) = (c + cir)re 



7 



(46) 



where ao still takes the value in Eq.([|5|), cq is determined by the normalization condition 
and c\ is a new variational parameter. By using variational method, we have 



E = 3.36202 



and 



v t = -0.01177 5u x = 0.00528, 
u 2 = -0.25348 5v 2 = 0.008526. 
Evidently, the accuracy of TWF is much more improved. 
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E. The other possible choice for the hypervirial operator 

One may choose different hypervirial operators G instead of what we employed above. 
Here, we give another example where the new hypervirial operator G' is in the form of 



With the Hamiltonian 
where 

we have the commutator 



G' = f(r)p 2 r . 



H = p 2 r + V eff (r) 



Veff(r) = ^ + V(r) 



in, g> } = m^M + mr)^ Pr -im pl - 2i m rf . (47) 

Let us set I = and V(r) = r for demonstration. Taking /(r) = r, we obtain a new first 
order hypervirial relation 

(-irpr) = {-ipD- (48) 

Similar to above operation, we define 

u[ = (-ipl) - (-irp r ) (49) 

and 

5u[ = j^r (50) 

(-« pf) 

as a criterion to judge the quality of a TWF. However, considering 

Pr = H — r, 

p r r + rp r = — i + 2rp r 
and Eq.(^), it is possible to give an alternative expression for the criterion, which reads 

v'l = E(-ip r ) + l-2(-irp r ). (51) 

Because for any eigenstate \R e i g ) of the Hamiltonian H we have 

H\R e ig) = E\R eig ), (52) 

v[ is identical to u" for the exact solution. For an approximate solution, these two expressions 
would give different results. However, it is not difficult to prove that for any function 
R(r) = u(r)/r which satisfies the normalization condition ( |A3| ) and the boundary condition 
we surely have 

(-iPr) = (53) 
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and 



(-irp r ) = -. 



(54) 



Consequently, v'[ vanishes trivially and cannot be used as a criterion. In general, v[ is not 
zero. Therefore, u[ and Si/[ can be applied to judge the quality of a TWF. We list the 
numerical result for the 15* state with the TWF parameters used in (|BT|) in Table 4. 

Table 4 



Variational results with the single-parameter TWF (Bl) for the 15 state in the linear potential case. 
The employed hypervirial operator is G". 



b 


1 


3 
2 


7 
4 


2 


v'l 


1 


0.16667 


0.018228 


-0.07559 


8u[ 


0.6667 


0.25000 


0.03517 


0.17810 



It is easy to see from Table 4 that 8u[ can indeed be used as a criterion to judge the 
quality of the TWF and gives a consistent result with those in Table 1, but it is less sensitive 
than £z/2 presented in Table 1. 



IV. CONCLUSION AND DISCUSSION 



There are a variety of approximation methods which can be employed to solve the non- 
relativistic Schrodinger equation where exact analytical solutions cannot be obtained. The 
variational method is one of the most frequently adopted approaches because of its apparent 
advantages. As well known, the obtained eigenenergies and wavefunctions are closely related 
to the form of the variational trial wavefunction (TWF) and the number of the variational 
parameters. What we concern is how close to the real values the achieved results can be. 
This is a crucial problem which should tell us the reliability of the results. In other words, we 
need a proper way to judge the quality of the TWF, especially the single-parameter TWF, 
because of its simplicity. 

To answer this question, we can set some criteria which determine the quality of the 
adopted TWF. In this work, we have carefully studied several calculable quantities for 
this purpose. Some of them were suggested in literatures and the rest of them are newly 
introduced by us. To obtain the idea of the criterion for judging the quality of TWF, we 
use the simple potential form V{r) = r, which has exact solutions for nS states, as an 
example. The numerical results indicate that all the criteria are consistent, namely, if one 
of them shows a sufficiently smaller deviation of the obtained result from the exact solution, 
the others confirm the conclusion. However, except the hypervirial relations, evaluating the 
quantities which stand for the criteria of the quality of TWF needs the exact solutions of the 
Schrodinger equation. Apparently, such solutions are generally not available. Only can the 
deviations from the hypervirial relations serve as a practical and powerful criterion for the 



17 



quality of TWF. Furthermore, in our special case with V(r) = r, we can deduce a conclusion 
that the hypervirial relations can be a good criterion and is supported by the other criteria 
gained in the special case where the exact solutions for nS states are available. This is 
also valid for various potential forms. To consolidate our confidence, we have studied the 
IP state with V(r) = r, where no exact solution exists. The results are qualitatively and 
quantitatively consistent with what we expected. 

Moreover, we use G = r N p r as the hypervirial operator which is of obvious physical 
significance. It should be specially noticed that in some special or accidental cases, the 
first (maybe the first a few) hypervirial relations precisely hold. However, it is not due to 
the correctness of TWF, but the TWF's fall in the special categories which are specified 
in Sec.III(D). In these cases, the first a few hypervirial relations cannot be used as criteria 
for the quality of TWF, instead one needs to invoke higher orders hypervirial relations or 
use other hypervirial operator, for example, G' in Sec.III(E). Thus, when one applies the 
hypervirial relations as criteria for the quality of TWF, he has to study the form of the 
adopted TWF to see if it falls into the categories listed in Sec.III(D). 

Our conclusion affirms that the hypervirial relations hold precisely for the real eigen- 
states, but approximately for TWF. They can be used to determine the quality of the 
adopted TWF. The deviation of the hypervirial relations from equality can serve as the 
most practical and powerful criterion for choosing an appropriate TWF in any potential 
cases. 
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APPENDIX A: THE ANALYTICAL SOLUTION OF THE SCHODINGER 
EQUATION FOR S-STATES IN THE LINEAR POTENTIAL CASE 

In the framework of the non-relativistic potential model, with a central potential V(r) = 
V(r), the eigenstate energy E and the corresponding radial wave function R(r) can be 
obtained by solving the Schrodinger equation 



vl Ki + i) 

2(i 2/ir 2 



H R{r) = ^- + K n ^J J + V(r) R(r) = E R(r), (Al) 



where H is the Hamiltonian, (i the reduced mass, p r the radial momentum operator and / 
the quantum number of the angular momentum. The radial wave function R(r) satisfies the 
normalization condition: 



oo 

2/„\ „2 



R\r) r 2 dr = 1 . (A2) 



Taking the linear potential V(r) = r/r^ normalized as ro = 1/GeV 2 (so we ignore tq in all 
later expressions), 2/x = 1 GeV and the reduced radial wave function u(r) = r R(r) which 
satisfies the normalization condition 



oo 

2/ 



u\t) dr = 1, (A3) 

o 

and the boundary condition 

u (0) = and «(oo) = 0. (A4) 

The equation for the S state has a simple form: 



2 



Rewriting the equation as 



d 

(-— + r)u(r)=Eu(r). (A5) 



(^ + (E - r))u(r) = 0, (A6) 



and comparing it with the Airy equation [12 



(^ - *) W(s) = 0, (A7) 

one easily obtains 

u nS {r) = N n Ai{r - E n ), (n = 1, 2, 3, • • • and r > 0) (A8) 

where Ai(^) is the Airy function, N n the normalization constant and E n the eigenenergy of 
the nS state in the linear potential case, which is just the value of the negative n-th node of 
the Airy function Ai(z). The numerical values of E n are E\ = 2.33810 • ■ ■, E 2 = 4.08794 ■ ■ -, 
E 3 = 5.52055 • • • etc., respectively. 

APPENDIX B: THE VARIATIONAL SOLUTION FOR THE S-STATES IN THE 

LINEAR POTENTIAL CASE 

For the ground state, i.e. the IS* state, in the linear potential case with 2/x = lGeV, the 
nomarlized trial wave function (TWF) with a single variational parameter is taken to be 

R twf (r) = N e~ ar \ (Bl) 

with the normalization constant 



= N e~ a r " 


yb(2a)U\ 


UttixdJ • 



N= v ' 2 . (B2) 

U7rr(§)J l ; 

In above equations, a denotes the variational parameter and h represents the model param- 
eter which specifies the type of the TWF. If one considers the linear potential only, the 
Hamiltonian for the nS state can be written as 
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d 2 2 d , nn . 

H = -—--- + r, (B3) 

ar z r ar 



Solving the Schodinger equation in terms of the TWF Eq.(Bl), one obtains the analytical 
results for a® and E. The optimal value of a is 



ao = 2( (i + 6)rm ) ' (B4) 



6 

and the corresponding minimum value of the energy reads 

E _ 3 r (i + 6)r(i)i*(jM rB ,, 

^--^f(|)l 2 ) ■ ( B5 > 

Because the behavior of the linear potential lies between the Coulomb potential and the 
harmonic oscillator potential, we select four kinds of TWF with b = 1 ( "hydrogen- like" 
wave function ) , b = |, b = % and b = 2 ( Gaussian wave function ), respectively. The 
resultant a^ and E min for various values of b are shown in Table 1. It is clear that when 
b = 7/4, one reaches the most accurate solution. 

Taking b = 7/4, one can construct a TWF for the 2S state. It reads 



R 2 t"f{r) = (co + c ir + c 2 r 2 )e- ao r * , (B6) 



where a = 0.348957, which is the optimal value determined for the IS state, cq can be given 
by the normalization condition, c\ should be determined by the orthogonal condition to the 
IS state and c 2 is the unique variational parameter which can be obtained by minimizing 
the corresponding average value of the Hamiltonian H, (R^JH\R^f). 
Similarly, for the 35* state, we can take TWF as 

Rtwf(r) = (co + c x r + c2r 2 + c 3 r 3 )e- a ° rl . (B7) 

The numerical values for the 2S and 3S states are also shown in Table 1. 
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